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1 . INTRODUCTION 


Two  separate  and  distinct  problems  in  the  modeling  of  turbu- 
lent flows  are  under  investigation  at  A.R.A.P.  under  AFOSR  support 
under  Contract  No.  F44620-76-C-0048 . 

1.  A method  for  extending  the  calculation  of  turbulent 
dynamics  by  second-order  closure  to  include  the  effect  of 
"structure"  as  well  as  scale. 

2.  The  development  of  a unique  closure  for  the  equations 
governing  molecular  mixedness  and  turbulent  chemistry  in  turbu- 
lent flows. 

Significant  progress  has  been  made  toward  the  solution  of 
problem  1 above  and  will  be  reported  in  Section  2 of  this  report. 
Although  the  results  are  very  new,  we  believe  that  a new  principle 
has  been  discovered  which  permits  rational  closures  for  turbulent 
mixing  and  turbulent  chemically  reacting  flows.  The  closures  are 
unique  and  satisfy  all  the  statistical  constraints  and  all  the 
limiting  known  conditions  that  we  have  been  able  to  check  the 
principle  against. 


2.  MODELING  TURBULENT  STRUCTURE 


In  current  second-order  closure  models  for  turbulent  flows, 
use  is  made  of  the  complete  equations  for  the  Reynolds  tensor 


uTuT  . These  models  contain  two  scales  associated  with  the 
turbulent  motion:  an  energy  containing  scale  A , and  a dissipa- 
tive scale  A . Since  the  large  scale  A is  related  to  the 
integral  scale  of  the  turbulence,  which  is  not  the  same  for  all 
velocities  and  all  directions,  our  equations  which  use  A in  the 
sense  of  an  average  scale,  do  not  allow  one  to  take  account  of 
turbulent  structure.  As  an  example  we  expect,  in  general,  that 
the  scale  associated  with  streamwise  velocity  fluctuations  (at 
two  points)  can  differ  significantly  from  the  scale  associated 
with  crosswise  fluctuations.  Further,  the  scale  depends  on  the 
orientation  of  the  two  points. 

This  information  is  fully  contained  in  the  tensor 


Rij(x,f)  - <u:(5  + I)  u;(S  - f)> 


In  particular: 

(i)  as  C ^ 0 , R^j  uTuT 

(ii)  The  integral  of  R^j  along  a direction  of  C yields 
the  integral  scale  corresponding  to  that  direction. 


We  propose  to  develop  and  analyze  a model  equation  for  the 

tensor  R^^  . It  follows  from  (i)  that  the  model  equation  for 

R^j  must  be  chosen  so  as  to  be  compatible  with  the  equation 

chosen  for  uTuT  . 

1 J 

Our  first  work  on  the  problem  of  structure  was  directed  at 
exhibiting  the  difference  in  the  integral  scales  associated  with 
streamwise  and  crosswise  components  of  velocity  fluctuations. 

We  thus  introduced  the  scale  tensor  A^j  which  measures  the 
scale  for  the  velocity  components  i and  j averaged  over  the 
orientation  of  the  two  points: 
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Rij  (x.D 
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We  further  define  the  mean  scale  by 

^ “ I ^kk 


Note  that  this  is  consistent  with  the  meaning  of  A as  it  is  now 
used  in  second-order-closure  calculations.  For  isotropic  turbu- 
lence we  can  express  A in  terms  of  G.I.  Taylor's  integral  scales 
as 


A 


L.  + 2L 

_f £ = 


. 666L£ 


and  it  is  surprising  to  find  that  a great  many  real  flows  are 
adequately  calculated  using  this  simple  result. 

To  obtain  the  governing  equation  for  the  scale  tensor,  we 
took  the  following  steps: 

1.  is  expanded  in  terms  of  its  moments  with  respect 
to  t . 

2.  The  expanded  equation  is  integrated  over  relative 
separation. 


The  resulting  equation  for  A^j  can  be  decomposed  without 
approximation  into  two  equations:  a scalar  equation  for  the  mean 
scale  A and  a (tracelen  tensor)  equation  for  the  deviator 


■ ''ij  - 3 \k«ij 


Since  the  equation  for  contains  the  source  q^A| 


3U. 


9U. 


9x.  + 9x.^ 


we  expect  that  the  deviator  will  reach  an  asymptotic  (super- 
equilibrium) value  on  scales  of  space  and  time  short  compared  to 
those  on  which  the  mean  scale  A will  evolve.  To  demonstrate 
this  new  result,  by  a special  example,  we  choose  the  simple  form 


This  new  result  and  the  results  which  follow  were  presented  at  a 
lecture  given  at  the  ONR  Martime  Meteorology  Workshop,  Boston, 
September  8-9,  1977,  by  Dr.  Coleman  duP.  Donaldson. 
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and  obtain  a pair  of  coupled  equations  for  the  mean  scale  A 
for  the  anisotropy  parameter  o . These  equations  are: 


and 


DA 


uTuT  9U, 


Dt  ■ ^ HT  + 

li^lqAf^  (q^A) 


V. 


q^  [ L 
1 9 


A 9x, 


qA  ’a!’’'' 


9x, 


+ 3'li-ir  + (3b  + 3V'  )q 


and 


The  parameters  and  b are  known  from  analysis  of  the 

Reynolds'  tensor,  Vg^  = -.03  from  observations  of  A in  grid 
turbulence,  p^  still  remains  to  be  determined  (we  have  given 
it  the  name  "ed^y  stripping  coefficient") . 

We  note  that  the  production  coefficient  in  the  A equation 

is  2^^  - ij  . Thus,  it  is  not  a constant  or  a flow  independent 

quantity  as  has  been  assumed  in  most  closure  calculations  made  to 

date.  For  parallel  flows  a popular  choice  for  this  nonconstant 

has  been  /_  x 

-2(1  - l)  = .35 

which  gives  o = 2.35  . 
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Since  parallel  flow  can  generally  be  estimated  by  super- 
equilibrium methods,  solving  the  o equation  under  the  assumption 
of  superequilibrium  yields 

a = 2.7 

which  is  in  surprisingly  good  agreement  with  the  value  currently 
popular.  We  note,  however,  that  a can  take  on  values  from  0 
to  3 , depending  on  the  type  of  mean  under  consideration.  The 
superequilibrium  values  for  the  streamwise  and  crosswise  scales, 
using  the  o just  obtained  are 


1.44 


A22  _ A33 

A A 


.78 


in  rough  agreement  with  observations, 


When  it  was  found  that  this  simple  example  showed  that  the 
"tensor"  character  of  the  scale  evolved  very  rapidly  compared  to 
the  mean  scale  Aj^j^/3  , it  suggested  that  the  same  result  would 
hold  for  the  two-point  correlation  tensor  , thus  permitting 

a far  more  detailed  analysis  of  the  turbulence  structure  than 
afforded  by  the  angular  averaged  A^j  . 

Following  this  line  of  reasoning,  we  separate  from  R. . the 
1 

scalar  quantity  R (which  determines  A)  and  obtain  the 


deviator 


^ij 


ij 


= R. 


4 R 6.  . 
3 aa  ij 


To  extend  current  computations  to  include  "structure"  to  a first 
approximation,  we  intend  to  determine  from  superequilibrium 

considerations  while  R will  be  calculated  dynamically.  To 
implement  the  latter  part  of  the  program,  we  only  need  a model  for 
two  terms  in  the  equation  for  R^^  that  vanish  when  t 0 and 
thus  do  not  appear  in  the  equation  for  the  Reynolds  tensor.  An 
attractively  simple  hypothesis  for  one  of  these  terms  (the 
cascading  term)  is  being  tested  against  grid  turbulence  data. 

This  work  is  continuing  and  is,  we  believe,  the  logical  way  to 
proceed  to  introduce  "structure"  as  well  as  scale  into  the 
calculation  of  turbulent  flows. 
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3.  MODELING  OF  MIXING 


During  the  reporting  period  we  have  further  developed  our 
models  for  scalar  variables.  We  have  obtained  a consistent 
second-order  closure  for  a binary  mixture  with  disparate  masses 
by  establishing  a box  model  which  is  completely  determined  by  the 
moments  calculated  by  the  second-order  closure  equations 
(a  , , p ‘ a ' p ' p ' ) , with  the  added  criterion  that  the  box 

model  should  correspond  to  the  most  unmixed  distribution  con- 
sistent with  those  moments.  The  criterion  is  motivated  by  three 
observations:  (i)  the  energy  containing  eddies,  which  are  of  large 
scale  are  now  thought  to  be  responsible  for  the  creation  of  new 
contact  surfaces  and,  if  this  is  so,  the  mixedness  will  be  less 
the  larger  the  eddy;  thus,  the  choice  of  minimum  mixing  consis- 
tent with  the  means  and  second-order  correlation  of  scalar 
qualities  that  are  being  computed  seems  appropriate;  (ii)  in  the 
presence  of  a fast  chemical  reaction,  it  is  essential  that  this 
particular  choice  be  made  if  a consistent  calculation  for  these 
fast  reaction  rates  is  to  be  achieved;  (iii)  this  particular 
choice  must  also  be  made  if  the  limiting  solutions  for  zero 
molecular  transport  are  to  be  recovered  from  the  closure  equations. 

Presently,  this  model  is  being  programmed  to  calculate  the 
experimental  results  being  obtained  at  CalTech  by  Roshko  and  his 
coworkers.  In  addition,  the  extension  of  these  ideas  to  more 
than  two  species  is  being  studied  theoretically. 
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